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Abstract 



Scattering on the PT-symmetric Coulomb potential is studied along a U- 
shaped trajectory circumventing the origin in the complex x plane from be- 
low. This trajectory reflects VT symmetry, sets the appropriate boundary 
conditions for bound states and also allows the restoration of the correct sign 
of the energy eigenvalues. Scattering states are composed from the two lin- 
early independent solutions valid for non-integer values of the 2L parameter, 
which would correspond to the angular momentum in the usual Hermitian 
setting. Transmission and reflection coefficients are written in closed analytic 
form and it is shown that similarly to other PT-symmetric scattering sys- 
tems the latter exhibit handedness effect. Bound-state energies are recovered 
from the poles of the transmission coefficients. 
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1 Introduction 



The Kepler-Coulomb problem has always played a special role in the for- 
mulation and the application of quantum mechanics. Besides being one of 
the textbook examples for exactly solvable problems, it also exhibits features 
that have always attracted the attention of mathematical physicists. Among 
these one can mention that the Coulomb potential possesses both discrete 
and continuous spectra, which can be associated with dynamical symmetry 
and Lie algebras describing them (see e.g. Ref. pQ for a review). Although in 
the description of realistic physical systems the three-dimensional Coulomb 
potential and the associated radial Schrodinger equation is used in most 
cases, much work has been done extending the Coulomb potential to other 
dimensions. Of these the one-dimensional Coulomb potential is the most 
notable, as the singularity at x = raises interesting questions both for the 
V{x) ~ —a;" 1 and the V(x) ~ — |a;| _1 potentials. (See e.g. Ref. [2] and refer- 
ences.) The discussion of this seemingly humble system requires techniques 
like the self-adjoint extension of the relevant differential operator [3]. 

Manifestly non-Hermitian versions of the Coulomb potential have also 
been studied in terms of PT-symmetric quantum mechanics. In this theory 
[I] Hamiltonians invariant with respect to the simultaneous V space and T 
time inversion were found to exhibit features characteristic for Hermitian 
systems, such as partly or fully real energy spectrum and the conservation 
of the norm. (See e.g. |5j for a recent review.) 

The first examples for PT-symmetric potentials were of the type V(x) = 
x 2 (ix) e , including the archetypal imaginary cubic potential for e = 1. After 
the first numerical results, the conjecture of the reality of the energy spectrum 
was proven analytically for such potentials [6] . An interesting aspect of these 
systems is that often they cannot be defined on the real x coordinate axis, 
rather their solutions are normalizable only along certain trajectories of the 
complex x plane. This was the case, for example, for the above potential 
with e > 2, when these trajectories had to fall into wedges lying in the lower 
half of the plane in e-dependent positions symmetrical with respect to the 
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imaginary axis. 

Later the PT-symmetric version of a number of exactly solvable poten- 
tials have been formulated, mainly along the real x axis or along a line 
parallel with it x — ic. (See e.g. Ref. [7J and references.) The importance of 
this imaginary shift was that singularities lying on the real x axis could be 
avoided, and at the same time, the energy spectrum remained independent 
of c. In this way real potentials defined on the positive half axis could be 
extended to negative x values too. As the result of this, solutions irregular 
in the Hermitian case became regular in the PT-symmetric version of the 
potential, and this led to a richer energy spectrum. 

The Coulomb potential was among the first exactly solvable potentials 
considered within the PT-symmetric setting. It was found, however, that it 
cannot be defined on the real X £1X1 S because the solutions were not regular 
for both x — > oo and x — > — oo [7J. In Ref. [8] a parabolic trajectory was 
proposed, which was inspired by the PT-symmetrized version of the well- 
known harmonic oscillator-Coulomb mapping. This study revealed that the 
spectrum of the PT-symmetric Coulomb potential includes a second set of 
discrete energy eigenvalues in addition to the one that is present also in the 
spectrum of the real Coulomb potential. But as a more interesting feature, 
the energy spectrum was inverted [8]. The interpretation of this unusual 
finding was given later in Ref. [9]. The transformation properties of the 
solutions under the PT operation including solutions both with real and 
with pairwise complex conjugate energy eigenvalues were also discussed [TO] . 
Another study of the PT-symmetric Coulomb potential was done in Ref. 
|llj : there the Coulomb potential was defined as V(x) ~ \x — ic| _1 . 

In the present work we extend, to the scattering scenario, the discussion of 
the "PT-symmetric Coulombic bound states as presented in refs. 0,19]. In this 
direction we intend to pursue two ideas. The first one reflects the existence 
of a number of publications [121 1131 EEU EES! EES] where the standard PT- 
symmetric version of the scattering problem has already been described and 
developed in application to a number of exactly solvable potentials. We feel 
inspired by the observation that in all of these works the scattering has only 
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been considered along the real x axis and/or along a trivially complexified, 
shifted straight-line contour x(s) = s — ic with real variable s and constant c. 
We shall change this perspective by employing an utterly nontrivial negative- 
mass generalization of the complex integration path x = xV £ Js) as already 
proposed, in the context of the stabilization of the Coulomb bound states, in 
ref. [9] (see also eq. (|5| below). In this setting we shall reveal the new role 
of the real parameter e which appears to bring a new degree of freedom in 
the phenomenological scattering theory, indeed. 

The second motivation of our present interest in the PT-symmetric Coulom- 
bic scattering along a U-shaped contour x(s) of ref. [9] is more physical since 
it reflects the unique possibility of the coexistence of discrete and scattering 
states in a single potential (in this respect cf., e.g., the review [TJ once more). 
We have to emphasize that in the PT-symmetric context such a feature has 
not yet been achieved even in the models of scattering along curved complex 
contours (cf., e.g., Ref. [T7] where the "tobogganic" integration path x = x(s) 
has been chosen even as extending, in principle, along several Riemann sheets 
of the scattering wavefunctions ip (x(s)) ). Thus, the present U-shaped choice 
of xV £ Js) will represent physics which varies with the "contour width" e. This 
seems to offer a scattering-scenario analogue of the variability of the bound- 
state spectra mediated, according to Refs. [U [TF], by the variability of our 
choice of the Stokes' "wedges" in the complex x plane. 

The structure of the paper is as follows. In Section [2] the general for- 
mulation of the problem is presented, together with the U-shaped trajectory 
along which the scattering problem is considered. Section [3] deals with the 
actual calculation of the transmission and reflection coefficients, while the 
results are summarized in Section |4j 



2 Definition of the problem 



Let us consider the Schrodinger equation 

h 2 



2m 



d 2 L{L + 1) 



dx 



x- 



ty(x) + V(x) $?(x) = EV(x) 



(1) 
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defined in the x variable, which runs along a trajectory of the complex x 
plane. Let us assume that this trajectory can be parametrized in terms of a 
real variable as x(s). In order to implement VT symmetry of this system we 
introduce 

V(x) = -, (2) 
x 

where Z and L(L + 1) are real. This latter condition is met if L is chosen 
real. The two linearly independent solution of ([I]) can be written in terms of 
confluent hypergeometric functions as 

qr^x) = Ci e- te s jL+1 lJ P 1 (l + L + iZ/(2k), 2L + 2, 2kx) (3) 

* 2 (x) = C 2 e- te a;- L 1 F 1 (-L + iZ/(2k), -2L, 2kx) , (4) 

where 2mE/h 2 = —k 2 . Note that ^ and Q represent the two linearly 
independent solutions only if 2L 7L. 

Equation Q with ([2]) differs from the usual radial Schrodinger equa- 
tion of the Coulomb problem in the complexified potential, and also in the 
trajectory it is defined on. In a systematic reformulation of real solvable po- 
tentials and the respective bound states to their PT-symmetric counterpart 
it was found [7] that this problem cannot be defined on the real x axis or its 
imaginary shifted version x — ic, because the boundary conditions cannot be 
implemented in both directions due to the exponential factor in (|3]) and Q. 
In an effort to determine the genuine PT-symmetric version of the Coulomb 
potential, the well-known Coulomb-harmonic oscillator mapping was used: 
this transformation was applied to the PT-symmetric harmonic oscillator 
defined on the imaginary shifted real axis x — ic j8]. The resulting trajectory 
was a parabola in the first and fourth quadrant, circumventing the origin 
from the left. In order to make it "PT-symmetric, i.e. left-right symmetric 
in the coordinate space, it had to be tilted to the first and second quadrant 
by the multiplication \x. In order to keep the kx quantity intact, k also had 
to be tilted in the opposite direction in the k wave number space as — ik. 
This resulted in the unusual finding that the energy spectrum was inverted, 
as can be seen from the relation 2mE/h 2 = —k 2 . Note that the two sets 
of discrete-energy solutions discussed in [8] are obtained from (|3| and Q 

6 



by substituting a non-positive integer in the first argument of the respective 
confluent hypergeometric functions, reducing them to the expected general- 
ized Laguerre polynomial form |18j . (In order to match the formulae of Ref. 
[H] with those obtained here the following substitutions have to be made: 
t -> x, A -> L + 1/2, e 2 = 1, while k 2 should be chosen Z/[2(n + L + 1)] 
and Z/[2(n — L)\ in the two cases, corresponding to the two possible values 
of the quasi-parity q = ±1.) 

As another possible trajectory, a U-shaped curve circumventing the origin, 
illustrated in the FigjTj was proposed in Ref. [9]. It is defined for a suitable 
e > as 



-i(s + f e) - e, se(— oo,— fe) 



a; (s 



i(s/e+3/27r) 



s G 



-£) 
2 ' 2 /' 



(5) 



i(s — |e)+e, s e (|e, oo). 
The asymptotic e— dependence of this curve has an immediate physical 




Figure 1: U-shaped curve in a complex x plane for e = 2 



meaning because it enables us to distinguish between non-equivalent alterna- 
tive asymptotes of the curves x(s) along which the non-equivalent asymptotic 
boundary conditions will be specified for our wave functions. Of course, in 
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contrast to the scattering wave functions which must be different on the left 
and right asymptotic branch of x(s), all of the curves of coordinates exhibit 
the same left-right symmetry x(—s) = —x*(s) in the complex x plane, which 
combines spatial reflection V with complex conjugation T that mimics time- 
reversal. It is seen that for large \s\ (in fact, for \s\ 3> sti/2) the solutions 
behave as exp(±ifc|s|). For real k = (—2mE) 1 / 2 /h (i.e., for m > and E < 
of Ref . jH] or for m < and E > in |Hj ) , this represents an oscillatory solu- 
tion. In parallel, for an imaginary k (i.e., for m > and E > or for m < 
and E < 0) it corresponds to exponentially decaying or growing solutions, 
depending on the sign of Im(/c). 

The next result of this analysis presented also in Ref. [§] is that at the 
negative mass, the U-shaped parametrization (JHl) opens the way towards the 
simultaneous description of the bound and scattering states. In this setting 
the role of the asymptotic physical coordinate is played by the real parameter 
s of course. 

The matching of logarithmic derivatives is usually used in models where 
the potential V is, at some point x(so) of the curve of (possibly, complexi- 
fied) coordinates, discontinuous. For analytic potentials (leading to analytic 
wave functions ip(x)), the situation is different since these functions are usu- 
ally well defined in all the points of some Riemann surface S. This means 
that in general, our analytic wave functions are multivalued functions which 
become single-valued, typically, on any selected Riemann sheet specified, 
say, as a subdomain D of a cut complex plane C. In such a scenario it 
is only necessary to match the logarithmic derivatives of our analytic wave 
functions ip(x(s)) during transition of the path x(s) between neighboring 
Riemann sheets of the Riemann surface (i.e., typically, between pairs of non- 
overlapping subdomains T>± of C). For this purpose it is usually sufficient to 
employ the analyticity of i/j(x) and to simplify the matching via a suitable 
deformation of the path x(s). Thus, most easily, one may analyse the tran- 
sition between T> + and T>_ just in an arbitrarily small vicinity of a branch 
point where functions ip{x) degenerate to their dominant parts with trivial 
analy t ic-cont inuat ion prop ert ies . 
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3 Scattering in the VT- symmetric Coulomb 
potential 

In what follows we shall make use of parametrization (|5]) to study scattering 
on the PT-symmetric Coulomb potential at negative mass. In Ref. [5] this 
unusual option has been explained as making the PT— symmetric Coulomb 
bound states stable. Here, we shall emphasize that such an option is also 
necessary for a consistent description of the scattering along the U-shaped 
complex contour. 

On a purely technical level we shall employ the natural analytic continuity 
of functions ^ and Q. In order to facilitate the implementation of this idea 
the auxiliary complex phase factors exp(27ri(L + 1)) and exp(27ri(— L)) will 
be introduced in the two solutions for Re (a;) > 0. Then, the asymptotic 
expansion of the solutions can be written as [19] 



l F l (a,b,z) ~ JgL( z -V-)- 2 F (a,l + a-6,-z~ 1 ) + 

m e * z a-b 2FQ ( b _ aA _ aiZ -i ) ) 

where lm(z) > 0, |Arg(^)| < n as \z\ — ► oo. 

Applying ^ to ^ and Q and employing the parametrization of x ^ 
the following asymptotic expansions are obtained for \s\ — > oo: 



where j '• = 1, 2. The logarithmic terms in the exponentials are characteristic 
of the Coulomb asymptotics and indicate that the Coulomb potential van- 
ishes slower than genuine short-range potentials exhibiting exponential tail, 
for example [20]. The coefficients are 

a l+ = Clla| ^ m Vie ^Jf^ (2t)) (8) 

01 . = aw-Mwn^ j^ ( 9) 

fc + = C l( 2 fc )--e^ e .,<3 I+ 3, e - Se - te _IX^)_ (10) 

= fliar^^ ^g, (ID 



d2+ 


= C 2 (2k) L e~ ikw£/2 e [7r2L e' 


*z k£ r(-2L) 

T(-L + iZ/(2k)) 


(12) 


a 2 - 


= C 2 (2k) L e ik7r£/2 e- [7vL e- 


nZ . T(-2L) 

T{-L-iZ/{2k)) 


(13) 


b 2+ 


= C 2 (2k) L e ik7r£/2 e- [7T3L e- 


-ke r(-2L) 

4fe e 

r(-L - iZ/(2fc)) 


(14) 


b 2 - 


= C 2 (2k) L e- [k ^ 2 e-^e 


- fa r(-2L) 

r(-L + iZ/(2fc)) 


(15) 



Note that as expected from the functional form of the solutions, the inter- 
change of the indices 1 <-> 2 corresponds to the interchange L <-» —L — 1. 
These coefficients are also related by the expressions 

6i+ = a 1 _e- 2fe e 2i7r ( x+1 ) &x_ = a 1+e - 2fce e- 2i ^ L+1 ) (16) 

6 2+ = a 2 _e- 2fc£ e 2i7r (- L ) 6 2 „ = a 2+ e~ 2k£ e 2inL (17) 

The asymptotic expansion of the general wave function $(x(s)) = aipi(x(s)) + 
f3ip 2 (x(s)) is then the following: 

$(s -> -oo) ~ (aa!_+/3a2_)e i(fcs -i ln( - 2fcs)) + («&!_+ /?6 2 _)e- i(fcs -i ln( ~ 2fes)) 

(18) 

$(s^ oo) ~ (aa 1+ + (3a 2+ )e [{ks+ ^ ln{2ks)) + (ab 1+ + (3b 2+ )b J+ e-^ ks+ ^ ln{2ks)) . 

(19) 

In the next step we may construct solutions that correspond to an incoming 
wave from one direction in order to evaluate the reflection and transmis- 
sion coefficients [16]. Using the coefficients above the following results are 
obtained after some manipulations with gamma, trigonometric and exponen- 
tial functions: 

T L ^ R (k) = ^-e-™ k£ e ^T(-L-iZ/(2k))T(L + l-iZ/(2k)) (20) 

R L ^ R (k) = T L ^«(A;)e- 2fc£ (-e-¥+2cos(27rL)) (21) 
T R ^ L (k) = -T L ^ R (k) (22) 
R R ^ L {k) = T R _ L {k)e 2k£ e=^ (23) 



It is seen that as expected, the bound-state energies emerge from the 



poles of Ti^ R {k) (20), i.e. when the arguments of the gamma functions are 



set to the — n non-positive integer. 
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Equation (22) is also in accordance with the results of scattering on other 
■pT-symmetric potentials in that the transmission coefficient does not show 
handedness effect, except for a trivial factor of —1, while the reflection coef- 
ficients clearly demonstrate handedness [TU [TH]. This means that similarly 
to other examples, waves arriving from the asymptotically absorptive side 
scatter differently from waves arriving from the asymptotically emissive one. 

To exemplify the results we present in Fig. [2] the potential as the func- 
tion of s for some parameter, together with the transmission and reflection 
coefficients in Fig. [3] for a series of parameters. As expected, the real and 
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(c) Z=1,L = 3.10, e = 0.005 



(d) Z=1,L= 3.01, e = 0.005 



Figure 2: The PT-symmetric Coulomb potential plotted as the function of 
s parameter for Z = 1, e = 0.005 and some values of L. 



the imaginary potential components are even and odd functions of s, further- 
more, the real component decays quicker than the imaginary one as it should 
based on Q. The role of the e parameter is essentially that of a length 
scale. L and Z set the scale of the real and imaginary potential component, 
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respectively. 
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(b) Z= 1,L = 3.53, e = 0.005 
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(c) Z = 1,L = 3.10, £ = 0.005 
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(d) Z= l,i = 3.01, 



0.005 



Figure 3: Transmission and reflection coefficients corresponding to the pa- 
rameters used in Fig. |2j 



It is seen that the large-fc behaviour of T^^k) is relatively smooth, while 
that of the reflection coefficients is dominated by the exp(±2fce) factors. An 
exception for this behaviour is seen only near integer values of L, which are 
excluded from the present discussion. The Z "charge" parameter does not 
influence the results in an essential way. 



4 Summary 

Scattering on the PT-symmetric Coulomb potential was discussed on a tra- 
jectory of the complex x plane. This trajectory was a U-shaped curve circum- 
venting the origin from below, guaranteeing the asymptotical regularity of 
bound states. From the topological point of view it is similar to the parabolic 
path obtained from the application of the harmonic oscillator-Coulomb map- 
ping of bound states to the PT-symmetric setting [S]. It is also reminiscent 
to the trajectories obtained by similar regularity arguments for the power- 
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type potentials appearing in the first publications of PT-symmetric quantum 
mechanics [I] : there the allowed wedges for the solutions were located in the 
lower half of the complex x plane, while here the allowed domain is the 
upper half, corresponding formally to the inverted power (ix) _1 expression 
of the Coulomb potential. In addition to the Coulombic potential term a 
centrifugal- like term L(L + l)x~ 2 was also considered. Here integer values of 
L were excluded, because in that case the linearly independent solutions had 
to be defined in another way. 

Parametrizing the trajectory as x(s) in terms of the real parameter s 6 
(—00, 00) allowed the expression of the asymptotic solutions in a form familiar 
from the discussion of the real Coulomb potential. When parametrized in 
terms of s, the real and the imaginary potential components vanished as 
s~ 2 and i.e. the imaginary potential component dominates the problem 
asymptotically. 

The transmission and reflection coefficients were determined, from the 
asymptotic solutions, with special attention to the continuity of the solutions 
in the complex plane. It was found that the transmission coefficients for waves 
arriving from the two directions differ only in a trivial factor of —1, while the 
reflection coefficients exhibit manifest handedness. This is similar to other 
examples for scattering in PT-symmetric potentials, which, however, were 
defined on the real x axis or its trivial shifted version x — ic. The e parameter 
appearing in the definition of the U-shaped curve plays a role similar to the c 
parameter applied in Refs. [T2J [33] to shift the trajectory off the real 
Although it changes the potential shape, it does not influence the energy 
spectrum. 

The results showed strong dependence on the e parameter that sets the 
distance of the U-shaped trajectory from the positive imaginary axis, while 
the Z "charge" parameter did not influence the results in an essential way. 
It essentially sets the relative weight of the imaginary and real potential 
components. Dependence on the L parameter was also significant in that 
the transmission and reflection coefficients showed rapid variations close to 
the forbidden integer values of L. 
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As another aspect similar to other scattering problems, both Hermitian 
and PT-symmetric, the two sets of bound-state energy eigenvalues described 
in Ref. [8] could be recovered from the poles of the transmission coefficients. 
Furthermore, the reparametrization of the problem in terms of the s variable 
also allowed to recover the energy spectrum with the correct sign: E < for 
bound states and E > for scattering states. 
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